Full tomography from compatible measurements 
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We put forward a reconstruction scheme prompted by the relation between a von Neumann measurement and 
the corresponding informationally complete measurement induced in a relevant reconstruction subspace. This 
method is specially suited for the full tomography of complex quantum systems, where the intricacies of the 
detection part of the experiment can be greatly reduced provided some prior information is available. In broader 
terms this shows the importance of this often-disregarded prior information in quantum theory. The proposed 
technique is illustrated with an experimental tomography of photonic vortices of moderate dimension. 
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Introduction. The quantum state is a mathematical object 
that encodes complete information about a system jit]: once 
it is known, the outcomes of any possible measurement can 
be predicted. Apart from fundamental reasons, acquiring the 
system state is invaluable for verifying and optimizing experi- 
mental setups. For instance, in some protocols of quantum key 
distribution, the knowledge of the entangled state distributed 
between the parties greatly limits the abiUtyof a third party to 
eavesdrop on the communication channel 12|]. 

The reconstruction of the unknown state from a suitable set 
of measurements is called quantum tomography jst]. Over the 
past years, this technique has evolved from the first theoret- 
ical 141] and experimental I5i concepts to a widely acknowl- 
edged and fairly standard method extensively used for both 
discrete ]0>I3] and continuous ^ variables. 

In this work, we focus on measurement strategies for the to- 
mographic reconstruction, leaving aside data post-processing 
issues. In practice, a sufficient number of independent obser- 
vations must be included in the set of measurements so that 
all physical aspects of the measured system are addressed. 
When dealing with complicated systems, such measurements 
may be difficult to implement in the laboratory due to various 
physical and technical limitations on the available controlled 
interactions between the system and the meter 

The goal of this Letter is to present a method of generat- 
ing a tomographically complete measurement set from a sim- 
ple von Neumann measurement that is readily implemented 
in the laboratory. Obviously, a von Neumann measurement is 
not complete, as all the measured projections are compatible 
and hence provide information only about the same aspects. 
However, as we shall show here, things are radically differ- 
ent when only a part of the full Hilbert space is of interest: In 
this subspace, even a simple von Neumann projection may be- 
come informationally complete. This should not be taken as 
an approximation, in the sense that some accuracy is traded for 
experimental feasibility. First of all, the energy of any system 
is always bounded, so one can restrict the attention to the sub- 
space spanned by low-energy states. Second, due to the finite 
resources, all quantum systems are de facto discrete and may 
be represented by a relatively small number of parameters. In 
that case, there is no necessity of sophisticated measurements 
that are informationally complete in the original large Hilbert 



space: since only a small subset is accessible, even much sim- 
pler observations are able to supply the information needed. 
This is the main idea behind the present contribution. 

Quantum tomography. Let us consider a density matrix g 
describing a c?-dimensional quantum system. A convenient 
representation of g can be obtained with the help of a trace- 
less Hermitian operator basis {Fi}, satisfying Tr(Fi) = and 
Tr(F,F,) = 5. 



Q = 



(1) 



where {a,;} are real numbers. The set {r^} coincides with 
the orthogonal generators of SU(o?), which is the associated 
symmetry algebra. 

In general, the measurements performed on the system 
are described by positive operator- valued measures (POVMs), 
which are a set of operators {Ilj} (with Hj > and J2j = 
1), such that each POVM element represents a single output 
channel of the measuring apparatus. The probability of detect- 
ing the jth output is given by a generalized projection postu- 
late = Trigllj). 

By decomposing the POVM elements in the same basis 
{Fj}, we get 



E 



(2) 



where {bj} are again known real numbers and C = {cji} is a 
real matrix. 

Informational completeness. A set of measurements will be 
called informationally complete if any quantum state g is un- 
ambiguously assigned to the corresponding theoretical prob- 
abilities pj. Since the projection postulate can be rewritten 

as 



(3) 



informational completeness requires the matrix C to have at 
least (P — 1 linearly independent rows. Numerically, this can 
be easily verified by calculating the rank of C, given by the 



2 




number of nonzero singular values. These are readily com- 
puted from the singular value decomposition of C. Thus, a 
set of measurements is informationally complete provided 



rank C> - 1 , 



(4) 



For example, a light mode can be treated as a harmonic os- 
cillator. The eigenstates of the rotated quadrature operators 
Q{9) = xcosO + psvL\9 comprise an informationally com- 
plete POVM. Naturally, only a finite set of projections can be 
done, so that a truncation of the original infinite-dimensional 
Hilbert space is necessary MM- In consequence, consider 
a von Neumann projection defined in the infinite-dimensional 
space Ti: X^fcLo ~ ^' where |fc) is an orthonormal ba- 

sis. Experimentally such measurements do not pose any dif- 
ficulty: all that has to be done is to determine the spectrum 
of a single observable. Nevertheless, this simple von Neu- 
mann measurement is not informationally complete in Ti, for 
all the observations are in this case mutually compatible and 
consequently no information about any of the existing com- 
plementary observables is available. 

Generating informationally complete measurements. As 
we will now show, this interpretation no longer holds when 
only a subspace 5 of Ti is considered. Let us specify S by 
introducing the projector Ps — X]f=o V''){^\ ' where |s) are 
eigenstates of Ps and S is the dimension. By projecting the 
original measurement on S, a POVM is induced in this sub- 
space, namely 



Y,nk = Y.Ps\k){k\Ps = t, 



(5) 



whose elements, in general, no longer commute [life, life'] 7^ 
0. Indeed, since the original commuting projections have dif- 
ferent overlaps with the subspace S, their mutual properties 
(commutators) are not preserved. In this way, an informa- 
tionally complete POVM may be generated. Obviously, this 
observation has many potential applications beyond tomog- 
raphy, although, due to strict space limitation only that topic 
will be discussed. 

The protocol we propose consists of the following steps: (i) 
A reconstruction subspace S is selected according to the par- 
ticular experiment, in such a way that all the relevant states are 
included, (ii) An experimentally feasible von Neumann pro- 
jection is chosen, (iii) The effective POVM induced in S, as 
given by Eq. (|5]l, is calculated and its informational complete- 
ness is checked with the help of condition If the induced 
POVM is informationally complete, the task is finished, oth- 
erwise the whole procedure is repeated with different choices 
of either the von Neumann projection or the reconstruction 
subspace or both. 

Before we proceed further, let us comment on the differ- 
ences between our protocol and the Naimark extension 11211 . 
which is another way of representing POVMs by projective 
measurements. This extension works by enlarging the Hilbert 
space with an ancilla, so the projective measurement acts on 
the product space of the system and ancilla. In our approach. 
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no. 1: Incompatibility (computed as the norm of commuta- 
tor) of the detections at two spatially separated pixels of a 
CCD camera in a truncated Hilbert space p = 0, . . . ,Pcutoff , 
i — — i?cutoff, ■ ■ ■ i^cutoff. Black (white) color means compatible 
(strongly incompatible), respectively. 



the possibility of representing a tomographical scheme by a 
projective measurement stems from the available prior infor- 
mation. In fact, the unpopulated states or unused range of 
variables play the role of ancilla here and, consequently, the 
measurement acts on a sum rather than a product space. 

Optical vortices. As a relevant example, we use our proto- 
col for the tomography of optical vortices. As the wave func- 
tion (or density matrix) in quantum theory, any transverse dis- 
tribution of complex amplitude (or coherence matrix) can be 
decomposed in a complete basis; the Laguerre-Gauss modes 
being a very convenient one 

hGl{x,y) = {x,y\£,p)^r\'\L\^\{2r^)e-^\^'t (6) 

where = + and = arctan(j//a;) are polar coordi- 
nates in the transverse plane and Lf, is a generalized Laguerre 



polynomial. It is well known 11 1311 that LGp beams exhibit he- 
licoidal wavefronts that induce a vortex structure and carry or- 
bital angular momentum of M per photon. Suppose a photon 
has been emitted into a superposition of modes, and we need 
to identify the resulting state. In general, this is an involved 
task fl4l ITsi [l6ll requiring the use of complicated optical de- 
vices. However, provided that only beams with bounded vor- 
ticities (i.e., values of |^|) are considered, as it is usually the 
case, our protocol can be employed and an informationally 
complete measurement can be generated from a very basic 
one, such as a single transverse intensity scan that is easy to 
record. In the language of quantum theory, this intensity scan 
is just I{x, y) (X Tr(gi|a;, y) {x, y\), where x and y denote now 
the coordinates of a given pixel of the position-sensitive de- 
tector Although detections in any pair of pixels are always 
compatible, in a subspace with bounded vorticities noncom- 
muting POVM elements can be induced. This is illustrated 
in Fig. [U which shows the noncommutativity (incompatibil- 
ity) corresponding to the positions {x, y) = (0, 0), and (0, 1) 



FIG. 2: Experimental setup of vortex tomography by means of com- 
patible observations. 

[in the same units of Eq. (|6]l]. Truncating the Hilbert space 
at smaller vorticities typically leads to stronger noncommu- 
tativity, although some nonmonotonicity is also observed as 
oscillations of gray shades appearing from the top-right to the 
bottom-left corner 

Experiment. To demonstrate the potential of the procedure, 
a full tomography of an optical vortex field from a single in- 
tensity scan has been performed in a controlled experiment. 
The experimental scheme is shown in Fig.|2] The beam gen- 
erated by a He-Ne laser is spatially filtered by a microscope 
objective and a pinhole. After the beam is expanded and col- 
limated by a lens, it impinges on an amplitude spatial light 
modulator (CRL Opto, 1024 x 768 pixels) displaying a holo- 
gram computed as an interference pattern of the required light 
and the inclined reference plane wave. 

Light behind the hologram consists of three diffraction or- 
ders (—1,0,+!), which can be separated and Fourier filtered 
by means of the 4/ optical system consisting of the lenses Li 
and L2, and an iris diaphragm. The undesired 0th and -1st or- 
ders are removed by an aperture placed at the back focal plane 
of the lens Li. This completes the preparation of a given state 
of Hght. 

Finally, a collimated beam with the required complex am- 
plitude profile is obtained at the back focal plane of the sec- 
ond Fourier lens L2, where a transverse intensity scan I{x, y) 
is acquired by a CCD camera. In the image plane, each pixel 
detection can be approximated by a projection on the position 
eigenstates \x,y){x,y\. As it has been shown above, while 
such detections are compatible in the full infinite-dimensional 
Hilbert space, an informationally complete POVM is induced 
in a subspace of truncated vorticities. 

In our experiment the superposition 

l^> = ^(K=l,P = 0) + |£ = 2,p = 0)) (7) 

was prepared by letting an amplitude spatial light modulator 
to display an interference pattern of the transverse amplitude 
(x, 2/1'!') and a reference plane wave, as mentioned above. Re- 
sults for this state are shown in Fig. [3] The ideal intensity 
distribution in the detection plane I{x,y) oc |(a;,y|^')p is 
shown in the left panel. This should be compared to the corre- 
sponding noisy recorded image shown in the middle panel. 
Finally, the right panel shows the best fit obtained with a 
maximum-likelihood algorithm IitIi in the subspace p = 
and £ = 0, . . . , 4. The reconstructed 5-dimensional density 



FIG. 3: Experimental tomography of optical vortex fields. From left: 
ideal intensity distribution, measured intensity distribution, and the 
corresponding best theoretical fit of measured data. 



matrix is shown in Fig. H) Notice that, due to experimen- 
tal imperfections (such as a discrete structure of the spatial 
light modulator, detection noise, etc.), the reconstructed state 
slightly differs from the ideal one (typical fidelities in our ex- 
periment are F « 96%). In view of the complexity of the 
system and the simplicity of the experiment, we consider this 
to be a very good result. 

Given the promising performance of the proposed scheme 
in this proof-of-principle experiment, the natural question is 
whether an experimentally feasible von Neumann measure- 
ment (such as a single intensity scan by a CCD camera with 
possibly very fine resolution) would furnish an information- 
ally complete measurement for any reconstruction subspace. 
To get some insights into this problem, we consider two dif- 
ferent scenarios related to the experiment above (see Fig. |5]). 
In the first case, only photons with nonnegative vorticities are 
considered: the full tomography from a single intensity scan 
is always possible. In the second case, both positive and neg- 
ative vorticities are allowed. Here a single intensity scan fails 
to provide complete information. It is easy to see why: since 
the intensity profiles of the Laguerre-Gauss modes LGp and 
LGp ^ are the same, perfect discrimination between states with 
positive and negative vorticities is not possible. Interestingly 
enough, some information about the negative part of the an- 
gular momentum spectrum is still available (see, e. g., the 
crosses in the plots for the same truncation ^cutoff), as it is also 
obvious from the fact that the phases exp{i£<f>) and exp(— 




FIG. 4: Real (on the left) and imaginary (on the right) elements of 
the reconstructed density matrix. 
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system when some prior information is available. This prior 
information does not only bring about a quantitative improve- 
ment of our knowledge, but may also make feasible a no-go 
task. Based on this observation, an efficient protocol was 
sketched providing the full characterization of complex sys- 
tems from simple measurements. This was demonstrated in an 
experiment with photonic vortices. In our opinion, this con- 
stitutes an improvement that will have a significant benefit in 
the number of different physical architectures where quantum 
information experiments are being performed. 
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FIG. 5: Informational completeness of measurements on vortex 
beams generated by a CCD camera with 11x11 pixels. The number 
of independent measurements 11^ generated from those 121 CCD de- 
tections are shown by circles for different truncations of the Hilbert 
space Ps- The number of independent measurements required for 
a complete tomography in the same reconstruction subspace is in- 
dicated by crosses. The reconstruction subspaces are truncated as 
follows. Upper panel: p = 0, ^ = 0, . . . ,£cutoff; bottom panel: 

P — 0, £ — —icMtoB, ■ ■ ■ , ^cutoff- 



in superpositions like LGq+LGq and LGg ^+LGJ can be dis- 
tinguished via interference with the other mode. This partial 
information is however not sufficient for the full characteriza- 
tion of this part of the reconstruction subspace. Provided one 
wants to keep the simple intensity detection, it is always pos- 
sible to use a fixed unitary transformation prior to detection 
to optimize the scheme. For instance, by increasing angular 
momentum of the measured beam by ^'cutoff (using, e. g., a 
charged fork-like hologram) the reconstruction subspace can 
be moved inside the nonnegative part of the angular momen- 
tum spectrum. This example nicely illustrates the role of prior 
information in experimental quantum tomography. 

Conclusions. We have shown that simple compatible ob- 
servations may provide full information about the measured 
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